Mon. Not. R. Astron. Soc. 000, [T]-?? (2009) Printed 25 August 2009 (MN I^TeX style file v2.2) 



The disc-jet coupling in the neutron star X-ray binary 
Aquila X-1 

: V. Tudose/'2'3*t R.P. Fender,'^'^ M. Linares,^ D. Maitra^ and M. van der Klis 

^Astronomical Institute 'Anton Pannekoek', University of Amsterdam, Kruislaan 4-03, 1098 S J Amsterdam, the Netherlands 
' ^Astronomical Institute of the Romanian Academy, Cutitul de Argint 5, RO-040557 Bucharest, Romania 

^Research Center for Atomic Physics and Astrophysics, Atomistilor 405, RO-077125 Bucharest, Romania 
, ^School of Physics and Astronomy, University of Southampton, Highfield, Southampton S017 IBJ 



Accepted 2009 August 24. Received 2009 August 24; in original form 2008 November 26 



ABSTRACT 

We study the accretion/ejection processes (i.e. disc/jet coupling) in the neutron 
star X-ray binary Aquila X-1 via a multi- wavelength approach. We use in the radio 
band the publicly available Very Large Array archive containing observations of the 
object between 1986-2005, in the X-ray band the archival Rossi X-ray Timing Explorer 
data (Proportional Counter Array and High Energy X-ray Timing Experiment) be- 
tween 1997-2008, and in optical (R band) observations with the Small and Moderate 
Aperture Research Telescope System recorded between 1998-2007. In the combined 
data set we find three outbursts for which quasi-simultaneous radio, optical (R band) 
and X-ray data exist and focus on them to some extent. We provide evidence that the 
disc/jet coupling in Aquila X-1 is similar to what has been observed in black hole X- 
ray binaries, at least from the point of view of the behaviour in the hardness-intensity 
diagrams (the hysteresis effect included), when the phenomenology of the jet is taken 
into account. Although based on a very small number of observations, a radio/X-ray 
correlation seems to exist for this system, with a slope of q;=0.40 ± 0.07 [Fradio oc F^), 
which is different than the slope of a=1.40 ± 0.25 found for another atoll source, 4U 
1728-34, but interestingly enough is relatively close to the values obtained for several 
black hole X-ray binaries. No significant correlation is found between the radio and 
optical (R band) emissions. We also report a significant drop in the radio flux from Aql 
X-1 above an X-ray flux of 5 x 10~^ erg cm^^ s^^. This behaviour, also reported 
in the neutron star X-ray binary 4U 1728-34, may be analogous to the suppression 
of radio emission in black hole X-ray binaries in bright, soft X-ray states. It suggests 
that from this point of view neutron star X-ray binaries can mimic the behaviour of 
black hole X-ray binaries in suppressing the jet in soft/disc-dominated X-ray states. 

Key virords: accretion, accretion discs - stars: individual: Aquila X-1 - ISM: jets and 
outflows - radiation mechanisms: non-thermal - X-rays: stars. 



1 INTRODUCTION 

The low-mass X-ray binary (LMXRB) Aquila X-1 (Aql X- 
1) is a recurrent soft X-ray transient which sho ws quasi- 
periodic outbursts about once a ye ar (Pricdhorskv fc Terrelll 
ll984l : lKitamoto et all [19931 : ISimonll2002f ). The compact ob- 
ject in the system is a neutron star as implied by the de- 
tection of type I X-ray bursts (e.g. iKovama et al.l Il98ll : 
ICzerny et al.lll987l : lYu et al.lll999l ). The companion is a main 
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seque nce star likely of the spectral type K7 (|Chevalier et al.l 
I1999I I. The X-ray timing and spectral properti es of Aql X - 
1 place it in the class of atoll sources ( Reig et al.l I2OO0I ) 
and recent observations of coherent pulsati ons in the per- 
sistent X-ray emission jCasella et al.l I2OO8I ) makes it the 
ninth accretion-powered millisecond X-ray pulsar (AMP) 
kno wn to date. The atolls show two main X-r ay states 
( Hasinger fc van der Kli^ll989l : Ivan der Klislll989l ), identifi- 
able in the colour-colour diagrams: a softer, "banana state" 
(BS), and a harder, "island state" (IS). In order to accom- 
modate later observat ions of harder X-ray spec tra a new 
state was defined (e.g. iPrins fc van der Klis|[l997l ): the "ex- 
treme island state" (EIS). From the point of view of the 
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X-ray spectral and timing properties, atoll sources (espe- 
cially in the EIS) share many properties with t he black hole 
XRB s (BHXRBs) in t he low-hard state (e.g. van der 
1 1994 lOlive et al.lll998l : iBerger fc van der K Us 1998,). 

A ql X-1 also exhib its X-ray burst oscillations 
IZhang et al.l Il998l ) and kHz quasi-periodic 
(QPOs ) both lower kHz QPQ (e 



(e-g- 

oscillations 



Qui et all Il998l: iMendez. van der Klis fc FordI [ioo 



Reig. van Straaten fc van der Klis 2004), and upper 

kHz QPQ (a single detection, iBarret. Boutelier fc Miller! 
l2008l ). 

The optical counterpart has in quiescence a V band 
magnitude of 21.60 and is only 0'.'48 away from a contam- 
inatin g star of V band magnitude 19.42 (jChevalier et al.l 
1 19991 ) thus complicating the optical studies. 

Besides the orbital period close to 19 h 

JChevaUer fc I lovaiskv 1991, 1998; IShahbaz et aD 1 19981 : 
IWelsh. Robinson fc Youndl2000i ). little is known about the 
other parameters of the system. The inclination of the orbit 
is po orly constrained, with values ranging between 36° and 
70° (Shahbaz, Casares fc Charles! [1997I : ICarcia et al] 1 19991 : 
IWels h. Robinson fc Young 2000i ). The distance to Aql X-1 
is not well kno wn, but different estimates place it in the 
range 4-6.5 kpc dChevalier et al ] |l999l : iRutledge et"alll200ll : 
IJonker fc Nelemansll20oi K 

Reports of radio obse rvations of Aql X-1 are extremely 
scarce in the literature: iHiellming. Han fc Roussel-Duprd 
l|l990l ): iRupen. Mioduszewski fc DhawanI (j2004l . l2005h . Tb 
date, only five radio detections have been documented, 
all obtained by the VLA during outbursts: four at 8.4 
GHz (two in 1990 August and two in 2004 May), and 
one at 4.9 GHz (in 2005 April). This is not surprising 
since the atoll sources are quite difficult to observe given 
their sub-mjy flux density levels at cm wa velengths even 
during outbursts (Fonder fc Kuulkerd I2OO1I ). Basically, al- 
though they represent the bulk of the XRBs population 
(jFendeJlioO^ ). only a handful of ato ll sources h ave been 
detected in radio (e.g. Hiellming fc H an 1995; M artf et al.l 
Il998l : iMigliari et al.l[2003l . l2004l : iMoore et al.ll200oi ). In the 
case of the atoll source 4U 1728-34 a correlation was 
found between the radio and X-ray fluxes in t he hard 
state (jMigliari et all l2003l : IMigliari fc Fendeill2006l ). resem- 
bling the relation established for BHXRBs in the low- 
hard state teorbel et al."2003':'Gallo. Fender fc Pool evll2003l : 
iGallo et al..,2006: Co rbel. Kording fc Kaarct 2008, but see 
also lXue fc Culf2007i ). but with a steeper gradient. 



2 OBSERVATIONS 



2.1 Radio 



We have analyzed all the public Very Large Array (VLA) 
radio data between 1986 February and 2005 December con- 
taining observations of Aql X-1 (close to 100 epochs). The 
majority of the observations were carried out at multiple 
frequencies, particularly at 4.9 and 8.4 GHz. Other frequen- 
cies, namely 1.5 and 15.0 GHz, were observed much less 
frequently. The bandwidth used was 100 MHz. All the runs 
were performed after reports of increased activity from Aql 
X-1 mainly at optical and X-ray wavelengths, and thus gen- 
erally trace the outburst history of the object. Given the 



target of opportunity nature of the observations, the ar- 
ray configurations were varying between runs, and the ef- 
fective observing time ranged roughly between 5 and 60 
min, with the median around 10-20 min. Primary cali- 
brators were either 3C286 or 3C48. As secondary calibra- 
tors we used, upon availability, J1950-I-081 or J1925-I-211 
for the vast majority of the epochs, J1922-I-155, J1939-100, 
J2007-h404, J1851+005, or J1820-254. The data were cali- 
br ated in AIPS04 using standard procedures as highlighted 
Diamoridlll 19951) . The imaging was done in DIFMAP 2.4e 



l|Shepherdlll997i r by progressively cleaning t he res idual map 



using the CLEAN algorithm of iHogbonj ([T97I). The uv- 
plane fitting was carried out with elliptical Gaussians within 
the same software. 

We detected Aql X-1 at more than 3ct levels on 7 epochs 
at 8.4 GHz (Table 1) and 4 epochs at 4.9 GHz (Table 2). 
In two occasions the detections were quasi-simultaneous at 
both frequencies. At least in the cases of some outbursts it 
is likely that the lack of detection in the radio band is partly 
due to the very short integration times which generated rela- 
tively high noise levels. In Tables 1 and 2 we also include the 
upper limits (3 a) of the flux densities for the observations 
pertaining to outbursts for which there are radio detections. 

When detected, the object was unresolved. All the flux 
densities reported here are determined via ut;-plane fitting 
and the corresponding errors are estimated based on the rms 
noise of the images and the differences between fiux densities 
measured in the uv-plaxie using different initial conditions. 
Two contaminating sources are present in the field of view, 
~40 and ~290 arcsec away from the position of Aql X-1, and 
affect the quality of the images in compact configurations at 
lower frequencies (i.e. 4.9 GHz). However, this should not 
significantly influence the flux measurements. Our values 
are in agreemen t , with in th e errors, with those report ed by 
iHieUming etHI (|l990l ) and iRupen et all l|2004 l2005h . The 
only exception is the data at 4.9 GHz from 2005 April (MJD 
53465.54) when we clearly detect the target while in the pre- 
vious work only a marginal detection is noted. This inconsis- 
tency comes probably from a misidentification of the source 
on this particular epoch. Generally speaking, the radio flux 
density levels of Aql X-1 are relatively low. At times, noise 
features are present in the radio maps and can be misidenti- 
fied as the target. This could have happened on this instance 
we mention. 

2.2 X-ray 

We used archival Rossi X-ray Timing Explorer (RXTE) Pro- 
portional Counter Array (PCA) and High Energy X-ray 
Timing Experiment (HEXTE) data of Aql X-1 available in 
NASA's High Energy Astrophysics Science Archive Research 
Center (HEASARC). 

We filtered out data taken when the spacecraft was 
pointing away from the source (offset > 0.02 degrees), too 
close to the Earth (elevation < 10 degrees) or during and up 
to 15 minutes after the South Atlantic Anomaly passages. 
All RXTE products were extracted using the standard HEA- 
SOFT tools (version 6.5). 

In the case of the PCA, the count rates were extracted 
using the Standard 2 data of all active PCUs (with 16 s time 
resolution) and ave raged over each observation. Following 
previous work (e.g. iReig et al.ir2004h we used the following 
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Table 1. The VLA 8.4 GHz radio detections and a selection of upper limits (3 a) together with the corresponding quasi-simultaneous 
observations in optical (R band) and X-ray (PCA/RXTE and HEXTE/RXTE) bands. 



Date 


MJD 


Flux density 


MJD 


Flux density 


MJD 


Flux 


MJD 


Gount rate X-ray 




radio 


radio 


optical 


optical 


PGA 


PGA 


HEXTE 


HEXTE state 




[day] 




[day] 


[MJy] 


[day] 


[10^'' erg cm^^ ^^^] 


[day] 


[c s-i] 



1990 August 


48111.23 


415 ± 38 






1990 August 


48129.23 


150 ± 21 






2002 March 


52334.62 


< 69 


52334.39 


372 ± 12 


2002 March 


52355.67 


179 ± 26 


52354.87 


302 ± 10 


2002 April 


52379.60 


< 68 






2002 May 


52399.58 


< 57 


52400.82 


16 ± 2 


2002 May 


52413.52 


< 63 


52414.29 


1 ± 2 


2004 March 


53073.71 


< 183 






2004 March 


53077.48 


< 147 






2004 March 


53085.56 


< 174 






2004 March 


53091.56 


< 150 


53092.38 


297 ± 10 


2004 May 


53144.50 


179 ± 22 


53144.36 


610 ± 19 


2004 May 


53151.42 


216 ± 19 


53151.35 


793 ± 24 


2004 June 


53162.37 


205 ± 19 


53162.38 


540 ± 17 


2004 June 


53170.35 


274 ± 22 


53169.40 


440 ± 14 


2004 June 


53176.29 


< 114 


53177.15 


121 ± 5 


2004 June 


53183.29 


< 126 


53182.27 


51 ± 3 


2004 July 


53187.23 


< 117 


53187.30 


15 ± 2 



52335.04 


5.57 


± 0.56 


52332.92 


1.1 ± 0.3 


BS 


52355.84 


1.40 


± 0.20 


52356.69 


2.3 ± 0.4 


IS 


53074.90 


1.26 


± 0.11 


53074.91 


11.3 ± 0.3 


EIS 


53078.84 


1.06 


± 0.11 


53078.84 


8.7 ± 0.3 


EIS 


53085.93 


1.06 


± 0.11 


53085.93 


9.2 ± 0.3 


EIS 


53088.90 


0.84 


± 0.21 






EIS 


53144.63 


1.44 


± 0.14 


53143.65 


16.1 ± 0.7 


EIS 


53151.71 


1.97 


± 0.13 


53151.71 


17.1 ± 0.4 


EIS 


53162.67 


2.41 


± 0.20 


53162.67 


17.3 ± 0.4 


EIS 


53170.40 


3.89 


± 0.51 


53172.37 


0.6 ± 0.5 


BS 


53176.92 


0.40 


± 0.04 


53176.92 


4.2 ± 0.4 


EIS 


53182.75 


< 


0.04 






EIS 



Table 2. The VLA 4.9 GHz radio detections and a selection of upper limits (3 a) together with the corresponding quasi-simultaneous 
observations in optical (R band) and X-ray (PGA/RXTE and HEXTE/RXTE) bands. 



Date MJD Flux density MJD Flux density MJD Flux MJD Gount rate X-ray 

radio radio optical optical PGA PGA HEXTE HEXTE state 

[day] [niy] [day] [^tjy] [day] [10"^ erg cm-^ s^^] [day] [c s^^] 



1990 August 


48111.23 


289 ± 24 
















2002 March 


52355.67 


154 ± 17 


52354.87 


302 ± 10 


52355.84 


1.40 ± 0.20 


52356.69 


2.3 ± 0.4 


IS 


2002 April 


52379.58 


< 105 
















2002 May 


52399.54 


< 93 


52400.82 


16 ± 2 












2004 May 


53151.42 


< 255 


53151.35 


793 ± 24 


53151.71 


1.97 ± 0.13 


53151.71 


17.1 ± 0.4 


EIS 


2004 June 


53162.37 


< 276 


53162.38 


540 ± 17 


53162.67 


2.41 ± 0.20 


53162.67 


17.3 ± 0.4 


EIS 


2004 June 


53170.35 


< 246 


53169.40 


440 ± 14 


53170.40 


3.89 ± 0.51 


53172.37 


0.6 ± 0.5 


BS 


2004 June 


53176.35 


< 174 


53177.15 


121 ± 5 


53176.92 


0.40 ± 0.04 


53176.92 


4.2 ± 0.4 


EIS 


2004 June 


53183.29 


< 282 


53182.27 


51 ± 3 


53182.75 


< 0.04 






EIS 


2005 March 


53446.71 


< 159 
















2005 April 


53465.54 


245 ± 29 


53465.34 


444 ± 14 


53465.66 


1.13 ± 0.11 


53465.66 


10.4 ± 0.4 


EIS 


2005 April 


53472.50 


317 ± 30 






53472.81 


2.70 ± 0.47 


53472.81 


19.5 ± 0.4 


EIS 


2005 May 


53494.52 


< 156 


53492.33 


228 ± 8 


53493.79 


0.84 ± 0.07 


53493.79 


6.7 ± 0.3 


EIS 


2005 May 


53506.48 


< 207 


53508.31 


107 ± 5 


53505.85 


0.09 ± 0.01 






EIS 



energy bands: A: 2.0-3.5 keV; B: 3.5-6.0 keV; C: 6.0-9.7 
keV; D: 9.7-16.0 keV. The colours and intensity were de- 
fined as: soft colour SC=B/A, hard colour HC=D/C, inten- 
sity I=A-|-B-|-C-|-D. We applied dead-time corrections and 
subtracted the background using pcabackest (version 3.6), 
the latest models and updated SAA history file^. Follow- 
ing the recommended procedural, we used the faint source 
model when the source was below ~40 c/s/PCU and the 
bright source model otherwise. 

The PCA gains changed several times during the life- 
time of RXTE, generating different "gain epochs" . On top of 



^ http:/ /heasarc. gsfc.nasa.gov/docs/xte/pca_news. html 

^ http: / /heasarc . gsfc . nasa.gov /docs/ xte / pca_news . html#:quick_table 



these sudden change s, a long term gain decay is also present 
(|jahoda et al.|[2006l ). In order to confidently compare count 
rates measured at different epochs, months to years apart, 
we normalized the count rates to those of the Crab nebula, 
which is known to be a steady X-ray source. The conversion 
from count rate to Crab was done using the value of the 
Crab rate closest in time and within the same PCA gai n 
epoch Ce.g. iKuulkers et al.lll993 : Ivan Straaten et al.ll2003l '). 
To obtain the unabsorbed flux in the 2-10 keV band we ex- 
tracted background and dead-time corrected PCA spectra 
using only PCU 2, as it has been the most active and best 
calibrated throughout the RXTE mission. We estimated 
the background spectra and created response matrices us- 
ing pcabackest (as explained above) and pcarsp, respectively 
(both part of the same HEASOFT v. 6.5 package), and 
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added a 1% systematic error to all channels and grouped 
them in order to have a minimum of 20 counts per energy 
bin. We then fitted the spectra within XSPEC f f .3.2ag us- 
ing the following model: wabs*(bbody+diskbb+powerlaw), 
and fixing the hydrogen column den sity to 0.5 x f 0^^ cm~'^ 
ijChurch fc Balucinska-Churchll200lh . In one case the system 
was in quiescence and we determined the corresponding up- 
per limit on the flux using an absorbed power-law to fit the 
unmodeled background (diffuse emission). 

In addition, we extracted 20-200 keV dead-time and 
background corrected count rates from the HEXTE stan- 
dard mode data, and averaged them within each observa- 
tion. We averaged both rocking directions in order to mea- 
sure the background spectrum during each observation. For 
the analysis in section 6 we used cluster B data (excluding 
detector 2, which was damaged in March 1996), while for 
the rest of the paper cluster A data (all detectors). This 
choice was determined by the periods of anomalous rock- 
ing of cluster A, which otherwise might have affected the 
long-term analysis done in section 6. The overall gain and 
response matrix of the HEXTE detectors are independent of 
time and hence no Crab nebula normalization was applied 
to the HEXTE data. 



2.3 Optical 

For the current work we used part of the data from the 
long-term monit o ring campaign on Aql X-1 reported in 
iMaitra fc BailvrJ l|2008h . Details of the data reduction and 
analysis can be found in section 2.1 of the aforementioned 
reference. 

The R-band observations were taken using the 
ANDICA^fl instrument on the 1.0 m and 1.3 m Small and 
Moderate Aperture Research Telescope System (SMARTS) 
at Cerro Tololo Inter- American Observatory (CTIO). The 
data reduction was done via the standard IRAF pipeline 
developed for ANDICAM. 

Around 83 percent of the quiescent R-ba nd flux is es- 
timat ed to come from the contaminating star (| Welsh et al.l 
I2OOOI ) and was thus subtracted from the total observed flux 
from all the observations. 

For magnitude to flux conversion it was as- 
sumed that R=0 ma gnitude co rresponds to 3064 Jy 
llCampins. Rieke fc Lebo fskv 19851: iBessell. Castelli fc Pled 
1 19981 ). In lChevalier et al.l (|1999l ') the colour excess of Aql X- 
1 was estimated as E(B-V)=0.5 ± 0.1. The R-band extinc- 
tion Ar was determined using the standard total to selective 
extinction coefficient Ay / E{B — V)=<i.l and the interstel- 
lar extinction law A r/Av=Q.7A?, (|Rieke fc Lebofskvl [l985l : 
ICardeUi et al]ll989l ). 



3 MULTI- WAVELENGTH LIGHT CURVES 

For the three radio detections in 1990 August no multi- 
wavelength information was available and therefore we do 
not discuss these data here. We only note that during the 
quasi-simultaneous detections at 4.9 and 8.4 GHz on MJD 
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Figure 1. Multi- wavelength light curves for the outburst of Aql 
X-1 from 2002 March and the variation of the hardness ratio 
between two PCA bands (9.7-16 keV and 6-9.7 keV), and the 
HEXTE (20-200 keV) and PCA (2-16 keV) bands. If not ap- 
pearing in the plots, the errors are smaller than the size of the 
points. The open symbols in the radio light curve denote 3 cr 
upper limits. 



48111.23 the spectrum between these frequencies was in- 
verted (i.e. Q ^ where F^, oc 



3.1 2002 March outburst 

During this outburst (Fig. 1) Aql X-1 was detected in radio 
quasi-simultaneously at 4.9 and 8.4 GHz and showed then 
an inverted spectrum. The optical and PCA data clearly 
indicate that the source was passing through an active state. 
It is also apparent that the beginning of the outburst was 
missed. The level of the HEXTE count rate is low during the 
^40 days of monitoring and tends to rise slightly towards 
the end of the data set. As indicated by the hardness ratio, 
most of the period the system was in a soft state which it 
left suddenly and transited quickly to a hard state around 
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Figure 2. Multi-wavelengtli liglit curves for the outburst of Aql 
X-1 from 2004 May-June and the variation of the hardness ra- 
tio between two PCA bands (9.7-16 lieV and 6-9.7 keV), and the 
HEXTE (20-200 keV) and PCA (2-16 keV) bands. If not appear- 
ing in the plots, the errors are smaller than the size of the points. 
The open symbols in the radio light curve denote 3 cr upper limits. 

MJD 52355. The radio detections are coincident with this 
X-ray state transition. 

3.2 2004 May-June outburst 

This is so far, relatively speaking, the most well covered out- 
burst of Aql X-1 in the radio band (Fig. 2). The source was 
detected on 4 occasions and only at 8.4 GHz, however the 
lack of detection at 4.9 GHz might be partly attributable to 
the high rms noise due to the short integration times. On 
the epoch with the highest 8.6 GHz flux density detected, 
the quasi-simultaneous upper limit at 4.9 GHz allows us to 
constrain the spectrum in this frequency range to an in- 
verted one. The optical Hght curve is well sampled and the 
PCA data reveal a major outburst preceded by a smaller 
flaring event at about MJD 53150. After an apparent peak, 
the HEXTE registered a drop in the count rate soon after 
the onset of the intensity increase in the PCA band. This 
has been obs erved before in Aql X-1 during two previous 
outbursts, bv lYu et al.l l|2003l ). The hardness ratio diagram 
indicates that during the active period the system made a 



Figure 3. Multi-wavelength light curves for the outburst of Aql 
X-1 from 2005 April and the variation of the hardness ratio be- 
tween two PCA bands (9.7-16 keV and 6-9.7 keV), and the 
HEXTE (20-200 keV) and PCA (2-16 keV) bands. If not ap- 
pearing in the plots, the errors are smaller than the size of the 
points. The open symbols in the radio light curve denote 3 a 
upper limits. 



series of X-ray state transitions, from hard to soft and back 
to hard. A similar behaviour is observed for the smaller flare, 
when the hardness ratio decreased by 10 percent, with the 
difference that in this case the system never reached the soft 
state, but only temporarily softened its spectrum. The in- 
creases in the radio flux density seem to trace the intensity 
enhancements in the PCA band. 



3.3 2005 April outburst 

For this outburst (Fig. 3), observations at 8.4 GHz were 
not available. Aql X-1 was detected twice at 4.9 GHz. The 
optical data is quite sparse and the onset of the burst was 
not covered well. The PCA and HEXTE intensity peaks are 
more or less coincident. The hardness ratio diagram reveals 
that the system never left the hard X-ray state during the 
active period. 
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Figure 4. HIDs for the outbursts of Aql X-1 from 2002 March (upper-left; X-ray data between MJD 52321. 1-MJD 52366.9), 2004 May- 
June (bottom-left; X-ray data between MJD 53053.8-MJD 53182.8) and 2005 April (upper-right; X-ray data between MJD 53464.0-MJD 
53505.8). The cumulative HID for the all three outbursts is represented in the bottom-right panel. The hardness is the ratio of the count 
rates in the PCA bands 9.7—16.0 keV and 6.0—9.7 keV. The points highlighted in each HID are shown in the corresponding radio light 
curve. The other points do not correspond to quasi-simultaneous radio and X-ray observations. 



4 HARDNESS-INTENSITY DIAGRAMS 



The hardness-intensity diagrams (HIDs) corresponding 
to the three outbursts are represented in Fig. 4. The hard- 
ness ratio is the hard colour defined in section 2.2. Note 
that the scales for all the HIDs are identical. The system is 
tracing the HIDs counter-clockwise. 

In 2002 March (Fig. 4, upper-left) the radio upper limit 
denoted "A" corresponds to a soft X-ray state of the sys- 
tem at a relatively high X-ray intensity. This situation is 
reminiscent of the similar behaviour observed in BHXRBs, 
when the radio emission is "quenched" above some X-ray 
flux level. We postpone a further discussion of the issue to 
section 6. The system stays in a soft X-ray state for a while 
and then makes a fast transition to a hard state. The radio 
detection ("B") caught the system during this transitional 
phase. 

In 2004 May-June (Fig. 4, bottom-left) the system 
traced a full cycle in the diagr am. It looks simila r to the 
HIDs observed in BHXRBs (e.g. iFender et aT]|2004l ) and re- 



cently in a dwarf nova l|Kording et al.|[2008l ). The first four 
radio non-detections ("A" to "D") as well as the first three 
detections ("E" to "G") correspond to a hard X-ray state. 
Then suddenly the system jumped to a soft X-ray state and 
the radio flux density increased significantly ("H"), by more 
than 30 percent (with the errors on the individual measure- 
ments of less than 10 percent). On atimescale comparable to 
the rise of the radio flux density, the system dropped bellow 
the radio detection limit, in a hard X-ray state ("I", "J"). 

In 2005 April (Fig. 4, upper-right) the system stayed 
constantly in a hard X-ray state. It ex perienced what is 
sometimes called a hard outburst (e.g. iMigliari fc Fenden 
I2OO6I ). Th is has been observed in Aql X-1 on a few occa- 
sions (e.g. [Rodriguez et al.ll2006l ). 

Fig. 4 clearly illustrates hysteresis in Aql X-1 (the X- 
ray luminosity at the transition from hard to soft state 
is different, i.e. higher, than the luminosity at the transi- 
tion from soft to hard state). This kind of beh aviour was 
previously noted for Aql X-1 (iMaccarone fc C oppi 200d; 
iMaitra fc Bailvnll2004l : [Gladstone. Done fc Gierrinski..2007n . 
A HID with all the PCA data for all the X-ray outbursts of 
Aql X-1 between 1997 and 2008 clearly shows this (Fig. 5). 
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Figure 5. The cumulative HID for all the outbursts of Aql X-1 
observed with the PCA/RXTE between 1997-2008. A hysteresis 
phenomenon is clearly at work in Aql X-1. The hardness is the 
ratio of the count rates in the PCA bands 9.7-16.0 keV and 6.0- 
9.7 ke V. Superimpo s ed (se e section 7 for details) is the "jet line" 
of e.g. iFender et all ||2004|') . 



as the EIS, the softer one as the BS and the transitional 
state in between as the IS. 

An almost complete cycle in the CCD was traced during 
the outburst in 2004 May-June (Fig. 6, bottom- left). Before 
the flare the system was in the EIS. As the count rate in the 
HEXTE band started to decline, the PCA intensity went 
up triggering a fast jump of the system to the BS. Aql X-1 
stayed in this X-ray state for about 10 days and then, again 
relatively suddenly, within a couple of days, returned to the 
EIS via the IS. This kind of behaviour seems to be typical 
for m any of the outbursts of Aql X-1 (see also iReig et aU 
I2OO4I who analyze in some detail different outbursts than we 
do here). Although some outbursts do not undergo a full 
cycle in the CCD, like the one on 2005 April (Fig. 6, upper- 
right), they do follow the same general trend in the CCD. 
The general pattern is the following: right before the onset of 
a major outburst in the PCA band the system is in the left 
side of the EIS, migrates towards the right side of the EIS 
and as the PCA intensity increases considerably it jumps 
suddenly to the BS. Here, moving towards its leftmost side 
(sometimes with excursions back to the right side), the peak 
of the PCA intensity is reached. Then, as the PCA intensity 
decreases the system jumps back to the EIS, passing quickly 
through the IS. 



Three major flares have been relatively well sampled and the 
corresponding hard to soft transitions in Fig. 5 are there- 
fore not observationally biased. Thi s phenomenon of hys- 



teresis has been seen in BHXRBs (e.g iRodriguez et al.|[20o3 : 
iRossi et al.ll2004 iBelloni et al.ll200a : lDunn et al.ll2008 l') sug- 



gesting a common origin for the state transitions in NSXRBs 
and BHXRBs. In Fig. 5 the hard to soft transitions might 
appear to be significantly quicker than the soft to hard ones. 
This is an artifact due to the adding of the data together 
as some minor outbursts also follow tracks in the region be- 
tween 0.01-0.10 Crab. Whether or not they trace a scaled 
down version of the large-scale diagram is not clear, but 
the superposition creates a region with a higher density of 
points, which could falsely be interpreted as symptoms of 
slower soft to hard transitions. 



5 COLOUR-COLOUR DIAGRAMS 



6 CORRELATIONS 



Although established only for a limited sample of objects, 
a correlation seems to exist, extending over many orders 
of magnitude, between the X-ray and radio emissions of 
BHXRBs in the low-hard X -ray state (|Corbel et aU |2003|. 
20081 : iGallo et ahlbOOl I2OO6I ). Recent evidence (|Xue fc Cuil 



2007 : [Galldl2007l ) suggests that this correlation is not 



sal. More observations are definitely needed in order to settle 
the issue since given the relatively sparse nature of the data 
(particularly in the radio band) a confident conclusion can- 
not be reached at the moment. In the case of NSXRBs the 
situation is no better due to the fact that most of them sys- 
tematically show weaker radio emission than BHXRBs (e.g. 
iFender fc KuulkersI [200lh . However, a correlation between 
the radio and X-ray emission has also been found for some of 
these objects in cert ain X-ray states, in particular the atoll 
sourc e 4U 1728-34 (|Mighari et al.ll2003l : iMigliari fc Fended 
l200d ). 



The colour-colour diagrams (CCDs) for the three out- 
burst studied here in detail are shown in Fig. 6. The system 
is tracing the CCDs clockwise. 

Already in Fig. 6 the CCD for the three outbursts hints 
that Aql X-1 seems to spend most of the time in two main X- 
ray states during the active periods: a harder state roughly 
with HC>0.9, and a softer state with HC<0.6 and a nar- 
rower interval of variation of SC. This is clearly confirmed 
by the CCD for all the outbursts observed with the PCA 
between 1997-2008 (Fig. 7). Besides the two main X-ray 
states, a transitional one is now evident, linking the two. 
This CCD is, not surpris i ngly, entirely consistent with the 
CCD buih bv lReig et aU with all the PCA data be- 

tween 1997-2002 (note that we averaged the data over each 
observation) . Following these authors and the nomenclature 
for atoll sources, we have identified the harder X-ray state 



For Aql X-1 we have found indications in the data that 
also for this object a correlation seems to exist between the 
radio and X-ray bands. Given the small number of obser- 
vations under consideration this should be taken as tenta- 
tive only, subject to confirmation or refutation when more 
data will be available (particularly in the radio band). Fig. 
8 shows the radio fiux density versus the PCA X-ray fiux 
for Aql X-1 (data from Tables 1 and 2). When the system 
was in the EIS, the 8.4 GHz data (3 points only) reveal that 
the radio and X-ray bands have a positive correlation with 
a correlation coefficient of 0.78. The slope of the correlation 
function is a=0.29 ± 0.23 {Fradio «; Fx)- Here and further 
on, the errors on the correlation function are standard errors 
of the regression fit. When considering only the 4.9 GHz data 
in the EIS (2 points only) the slope becomes 0.30. Taking 
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Figure 6. CCDs for the outbursts of Aql X-1 from 2002 March (upper-left; X-ray data between MJD 52321. 1-MJD 52366.9), 2004 May- 
June (bottom-left; X-ray data between MJD 53053.8-MJD 53182.8) and 2005 April (upper-right; X-ray data between MJD 53464.0-MJD 
53505.8). The cumulative CCD for the all three outbursts is represented in the bottom-right. The soft colour is the ratio between the 
count rates in the 3.5-6.0 keV and 2.0—3.5 keV bands and the hard colour is the ratio between the count rates in the 9.7-16.0 keV and 
6.0-9.7 keV bands. 



into account all the radio detections in the EIS, independent 
of the frequency (i.e. 5 points), the data show no correlation 
between the radio and X-ray bands (correlation coefficient 
0.38). One possible explanation for this, assuming the cor- 
relation at 8.4 GHz is real, is that the spectrum between 
4.9 and 8.4 GHz was not flat during these observations. In 
BHXRBs for instance, in general in the low- hard state th e 
radio spectrum is flat or optically thick (e.g. lFendeij[200lh . 
however in a few occasions it was observed that the spectrum 
changes rapid ly and tends to bec ome optically thin prior to 
a r adio flare (iFender et al.l 120041 '). Perhaps something simi- 
lar happens in Aql X-1. Unfortunately, in our data set the 
only constraints to the radio spectrum comes from the ob- 
servations of the outbursts from 2002 March (Fig. 1) and 
2004 May-June (Fig. 3), when the system was not in the 
EIS. When all the 8.4 GHz data are taken into account, in- 
dependent of the X-ray states of the system (i.e. 5 points), 
the slope of the correlation function is q=0.40 ± 0.07 with 
a correlation coefficient of 0.96 (dashed line in Fig. 8). A 
similar analysis, for all the 4.9 GHz data (3 points) reveals 
an at best marginal correlation (correlation coefficient 0.60) 
with a slope q:=0.48 ± 0.64. 



The radio upper limit denoted "Ql" in Fig. 8 is notable. 
It has a relatively high PGA flux and represents a deep radio 
upper limit. The system was in the BS at the time (point 
"A" in Fig. 4, top left). It looks as if the radio emission was 
"quenched". A very similar situation is observed in some 
BHXRBs in the high-soft state, where the radio emissio n 
drops significantly above some X-ray flux l|Gallo et al.ll2003l ). 

4U 1728-34 is the atoll source with the most numer- 
ous quasi-simultaneous detections in the radio and X-ray 
band (i.e. 10 up to now). We have used these data (VLA 
at 8.5 GHz with 100 MHz bandwidth and PGA/RXTE in 
the energy band 2-10 keV) re ported in lMigliari et all (|2003l ) 
and lMigliari fc Fender! (|2006l ) to compare them against the 
Aql X-1 data from the point of view of the radio/X- 
ray correlations (Fig. 9). For Aql X-1, in calculating the 
luminosity we assumed the distance of 5.2 kpc used by 
iMigliari fc Fendeil (|2006l ). A flat radio spect rum extending 
up to 8.5 GHz was also assumed. Following IMigliari et al. 
(2003) and Migliari & Fender (2006,), using all the data for 
4U 1728-34 (except for the point denoted "Q2") we found 
a positive correlation between the radio and X-ray bands, 
with the correlation coefficient 0.90 and the slope of the cor- 
relation function 0=1.40 ± 0.25 {LraMo oc I/x), in perfect 
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Figure 7. The cumulative CCD for all the outbursts of Aql X-1 
observed with the PCA/RXTE between 1997-2008. While in the 
active state, the system spends most of the time in the "extreme 
island state" (EIS), with quick transitions through the "island 
state" (IS) and short excursions in the "banana state" (BS). The 
soft colour is the ratio between the count rates in the 3.5-6.0 keV 
and 2.0-3.5 keV bands and the hard colour is the ratio between 
the count rates in the 9.7—16.0 keV and 6.0-9.7 keV bands. 



Figure 8. The radio/X-ray correlation in Aql X-1. The filled 
circles and squares correspond to the radio detections at 8.4 GHz 
and 4.9 GHz respectively. The empty symbols are the radio upper 
limits (3 o"). The dashed line describes the correlation function 
(correlation coefficient 0.96) when taking into account all the 8.4 
GHz detections independent of the X-ray state of the system, and 
has a slope a=0.40 it 0.07 (iVadio F'^). The radio upper limit 
denoted "Ql" is discussed in section 6. 



agreement with previous authors. If moreover we consider 
only the data corresponding to the hard state (i.e. 7 points; 
see Mi gliari fc Fonder 2006 for an identification of the X- 
ray states) the correlation coefficient becomes 0.96 and the 
slope 1.42 ± 0.20. It seems therefore that the correlation 
function for Aql X-1 is flatter than the corresponding cor- 
relation function for 4U 1728-34. The normalization of the 
trends however, appears to be similar. This result seems to 
be significant, although given the very sparse data available 
at the moment in the radio band it is not possible to draw 
a solid conclusion. More radio observations are stringently 
needed in order to further constrain the correlation between 
the radio and X-ray ba nds for these sources . 

As noted before bv lMigliari et all (120031 ) the point "Q2" 
is interesting. It's the highest PCA X-ray flux at which the 
radio counterpart was detected, and the radio emission is 
relatively weak. Moreover, the system was perhaps in a soft 
X-ray state, but this is not clear given the slightly contra- 
dictory spectral and timing properties (see the discussion in 
iMigliari et alll2003l ). Thus, it might be that the radio emis- 
sion was "quenched" during this observation. The radio de- 
tection "Q2" of 4U 1728-34 and the radio upper limit "Ql" 
of Aql X-1 therefore suggest that NSXRBs might show a be- 
haviour similar to BHXRBs with respect to the suppression 
of the radio emission above some X-ray flux levels. 

No significant correlation was found between the radio 
and optical bands. The correlation coefficient for the 3 ob- 
servations in 2004 May-June for which quasi-simultaneous 
data exists at radio (8.4 GHz) and optical wavelengths (see 
Table 1) in the EIS is 0.52. Adding also the observations in 
which the system was in the IS and the BS does not improve 
the correlation coeffic i ent (0 .01). 

iMaitra fc Bailvnl (|2008l ) found correlations between the 
optical (R band) and ASM X-ray bands. The slope of the 
correlation function seems to be steeper when the system 



is in the hard state, however outliers are present and it's 
not entirely clear what role, if any, the observational bias 
is playing. In turn, we were able to confirm this behaviour 
using PCA instead of ASM X-ray data (Fig. 10, top). AU 
the multi-wavelength observations presented in Fig. 10 were 
quasi-simultaneous (separated by less than 1 day). The sys- 
tem in the EIS follows a different optical/X-ray correlation 
than in the BS. 

Aql X-1 also shows distinct behaviours in the EIS and 
the BS X-ray states when the optical (R band) and HEXTE 
bands are compared (Fig. 10, middle). In the BS the HEXTE 
X-ray count rate is systematically low, independent of the 
optical flux. This is contrary to the trend observed in the 
PCA band where it is in the EIS that the system shows low 
X-ray intensities at very different optical flux levels (Fig. 10, 
top). 

The interplay between the X-ray emission in the PCA 
and HEXTE bands can be seen in Fig. 10, bottom. Please 
note that in relation to Fig. 10 the term "correlation" was 
used loosely, more like in the sense of "trend". Formally 
speaking, due to the scatter in the data, the correlation co- 
efficients associated to Fig. 10 are not significant, with two 
exceptions: 0.93 for the optical/PCA X-ray bands in the BS 
(Fig. 10, top) and 0.95 for the HEXTE/PCA X-ray bands in 
the EIS (Fig. 10, bottom). Since only data during outbursts 
are available, and we only plotted the quasi-simultaneous ob- 
servations (thus only part of the total amount of data was 
selected), the analysis might be biased, in particular in the 
top and middle panels of Fig. 10. More quasi-simultaneous 
multi-wavelength observations are needed in order to con- 
strain the behaviour of Aql X-1. 



10 V. Tudose et al. 




T — ' — I — ' — I — ' — I — ' — r 



T — ' — r 




». 750 - 



O O no 



O 



O 

Q,0 oo o 




J 1 I 1 I 1 I 1 L 



• EIS 

O BS 
IS 



O.I 0.2 0.3 0.4 0.5 0.6 0.7 0.8 

PCA intensity [Crab] 



Figure 9. The radio (8.5 GHz)/X-ray (2-10 keV) correlations in 
Aql X-1 and 4U 1728-34. The filled circles and squares correspond 
to the radio detections of Aql X-1 and 4U 1728-34 respectively. 
The empty circles are the radio upper limits (3 a) for Aql X-1. 
The dotted line is the correlation function (correlation coefficient 
0.90) for 4U 1728-34 (independent of the X-ray state) and has a 
slope a=1.40 ± 0.25 (L^ 

adio ^ ^x)' The dashed line describes 
the correlation function (correlation coefficient 0.96) for Aql X-1 
when taking into account all the 8.4 GHz detections independent 
of the X-ray state of the system and has a slope o=0.40 ± 0.07. 
The data points denoted "Ql" and "Q2" are discussed in section 
6. 
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7 DISCUSSION AND CONCLUSIONS 

Given the scarcity of the radio observations of Aql X-1 in the 
literature, we have searched the pubhc VLA archive in the 
period 1986-2005 and found previously unreported detec- 
tions thus bringin g up the number of secure detections (more 
than 3 a) from 5 (iHiellming et al.lll990l : iRupen et aLll2004 
l2005t l to 11. To attempt a multi- wavelength characterization 
of the system we used archival RXTE X-ray data (PCA and 
HEXTE) between 1997-2008 and the op tical observations 
(R ba nd) between 1998-2007 reported in iMaitra fc Bailvnl 
l|2008h . 

In the combined data set on Aql X-1 we found 3 out- 
bursts for which quasi-simultaneous radio, optical and X-ray 
data exist (Figs. 1, 2 and 3). These light curves look quite 
different from each other, at all the wavelengths, showing 
that the source is revealing a rich phenomenology. Judging 
by the hardness ratio, the three outbursts under discussion 
can be considered as corresponding to distinct classes of out- 
bursts: in 2002 March the radio flare can be associated to 
a soft to hard transition; in 2004 May-June to a "classi- 
cal" hard to soft transitio n; and in 2005 April to a "hard 
flare". In the optical band lMaitra fc Bailvnl l|2008l) pointed 
out that the light curves of most of the outbursts of Aql 
X-1 do not belong to the pure fast rise exponential decay 
(FRED) topology and showed that FREDs occupy a partic- 
ular region in the optical/X-ray plane that is distinct from 
the non-FRED outbursts. 
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Figure 10. Trends in the behaviour of Aql X-1 in different X- 
ray states. Filled black circles correspond to the "extreme island 
state" (EIS), open circles to the "banana state" (BS) and filled 
gray circles to the "island state" (IS). All the observations in 
different bands are quasi-simultaneous (separated by less than 1 
day). The used energy bands are 2-16 keV for the PCA/RXTE 
and 20-200 keV for the HEXTE/RXTE. The optical band is the 
R band. 
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But in spite of the diverse morphologies of the light 
curves, the outbursts of Aql X-1 tend to converge to a com- 
mon behaviour when the HIDs and CCDs are investigated. 
The three outbursts presented in detail here, as well as most 
of the other events, track generally a similar path: counter- 
clockwise in the HIDs and clockwise in the CCDs. There 
might be nonetheless exceptions. In one case during a minor 
outburst in the soft state the system tracked partially the 
HID (Fig. 5, the left side of the diagram) clockwise. Dur- 
ing two or three minor outbursts the system was transiting 
between the hard and soft states and back again within a 
very small range in the X-ray intensity (Fig. 5, the bottom 
side of the diagram, between 0.01-0.10 Crab) and it's hard 
to say with certainty whether it mimicked a scaled down 
version of the diagram or was doing something completely 
different. The fact remains that the HIDs of Aql X-1 exhibit 
similar characteristics to the HIDs of some BHXRBs (e.g . 
Rodriguez et al."2003': ' Rossi et al.ll2004 iBelloni et aLll2005l : 



an open issue (e.g. iDone fc Gierliiiskill2003l: [ Zdziarski et al.l 



Dunn et al.l 2008. : .Fender et al. " I2OO9I ') . For these objects a 



phenomenological model was developed that accounts for 
the connection between the radio and X-ray em ission via 
the presenc e or absence of jets in the system (1 Fender et al.l 
I2OO4I . I2OO9I ). It goes as follows, starting from the lower right 
corner of the HID in Fig. 5. In the hard X-ray state a steady, 
optically thick jet is produced. During a typical outburst the 
X-ray intensity increases (moving up in the HID) and after 
a while the X-ray spectrum begins to soften (moving to the 
left in the HID). At the "jet line" the system enters the 
soft X-ray state when no jet is produced. The switching off 
of the jet is associated with relativistic ejections of matter, 
generally optically thin. The X-ray intensity decreases (mov- 
ing down in the HID) and soon the X-ray spectrum starts to 
harden (moving right in the HID). At the crossing of the "jet 
line" the system switches on the steady jet and finishes the 
cycle back in the lower right corner of the HID. This kind of 
behaviour is roughly consistent with what is observed dur- 
ing the three outbursts of Aql X-1 (Fig. 4), but given the 
small number of radio detections it is practically impossible 
to test the model in more depth. The radio detection in a 
soft X-ray state with a likely inverted spectrum during the 
outburst from 2004 May-June (point "H" in Fig. 4, bottom 
left) does not necessarily contradict the idea of a turned 
off jet in the soft X-ray state since ad-hoc explanations are 
readily available, for instance by invoking a blob of matter 
ejected during the hard to soft tr ansition at the "jet line", 
which is expanding adiabatically Ijvan der Laanlll966l ). 

The presence of hysteresis in Aql X-1 has alrea dy been 
estab l ished (e.g. Macca ronc & Coppi 2003; Maitra fc Bailvnl 
I2OO4I : [Gladstone et al.il2007f ) and confirmed in the present 
work (Fig. 5). Another atoll source, 4U 1608-522, shows a 
behaviour resembling that of Aql X-1 (Linares & van der 
Klis, in preparation). The similarity with the hysteresis phe- 
nomenon observed in BHXRBs (e.g. [R odriguez ct al. 200 3; 
iRossi et all |2004 iBeUoni et all l2005l : iDunn et all \2Qm ) 
seems to suggest a common mechanism. What might be 
its nature is still a matter of debate. This is connected 
to the fact that at the moment, no self-consistent physical 
model of the X-ray states exists, although a few promis- 
ing candidates are available (e.g. Remillard fc McClintockl 



l2004riMever-Hofmeister. Liu fc MeverlbOOSl ). 

Evidence for radio/X-ray correlations were found in Aql 
X-1 (Fig. 8). Taking into account the radio detections at 8.4 
GHz, the slope of the correlation function is q:=0.40 ± 0.07 
{Fradio oc Fx) with a correlation coefficient of 0.96. This 
is significantly different than the rad io /X-ray correlation 
for another atoll sourc e, 4U 1728-34 (Migliari et all l2003l : 
iMigliari fc Fenderll2006l ). for which the data indicate a slope 
of 1.40 ± 0.25 with the correlation coefficient 0.90, although 
the normalizations of the correlation functions are similar for 
both sources. However, these results are based on a very lim- 
ited sample of data points, 5 in the case of Aql X-1 and 9 for 
4U 1728-34, and therefore need further testing with future 
observations. Nevertheless, it's interesting to note that the 
slope of the correlations between the radio (8.5 GHz) and X- 
ray bands (2-10 keV) observed in some BHXRBs in the hard 
state are relatively close t o what we obtained for Aql X-1: 
0.51 ± 0.06 fo r V404 Cyg jCorbel et al.ll2008l ). 0.58 ± 0.16 
for A0620-00 (|Gallo et al.ll2006l 'l. But~the normalizations of 
the correlation functions are higher for BHXRBs. This is re- 
lated to the fact that the ratio between the radio and X-ray 
luminosity for BHXRBs is between one and two orders of 
magnitude higher than the corresponding ratio for N SXRBs 
iFender fc KuulkersI I2OOII : IMigliari fc Fendeij I2OO6I ). This 
can be explained for instance by a more efficient (i.e. "ra- 
dio loud") jet production m echanism in BHXRBs than in 
NSXRBs (jFender et al.ll2003l ) which might come about due 
to the fact that the geometrically thick disks are more ef- 
ficien t at producing outflows (e.g. iMeier. Koide fc Uchid3 
I2OOII ) and the disks of NSXRBs are truncated at larger radii 
than those of BHXRBs. Or it might be that in the case 
of BHXRBs a large fraction of the extracted gravitational 
energy is advected across the event horizon, the end prod- 
uct of the process b eing more "X-ray quiet" systems (e.g. 
iNaravan et al.iri997l ). From a more general point of view, 
it was shown that a scaling relation between the X-ray (2- 
10 keV) and radio (5 GHz) bands with a slope of 0.60 ± 
0.11 holds for a fairly large sample of BHXRBs and ac- 
tive galactic nuclei (AGN) wh en the mass of the compact 
objects is taken into account (IMerloni. Heinz fc di Matted 



I2OO6I : iDone. Gierliiiski fc Kubotall2007l for reviews). There- 
fore the trigger of the transitions between these states, which 
is likely the mechanism responsible for the hysteresis, is still 



I2OO3I : iFalcke. Kording fc Markofj 120041 '). Given also the 
cent advancements (e.g. the relation between the mass of 
the black hole, the accretion rate and the bre ak frequency 
in X-ray power spectra for AGNs and XRBs ^McHardv et al.l 
or the similarity between th e HIDs of XRBs and the 
HID of a dwarf nova in outburst [Kording et al. I l2008l ]). it 
seems that evidence is mounting at a fast pace that the ac- 
cretion/ejection process works in a similar way in a broad 
class of objects: stellar mass black holes in XRBs, super- 
massive black holes in AGNs, and perhaps neutron stars in 
XRBs and white dwarfs in cataclysmic variables. 

Aql X-1 shows evidence for the "quenching" of the ra- 
dio emission in the soft X-ray state, a bove some X-ray flux 
(the radio upper limit "Ql" in Fig. 9). IMigliari et al.1 (|2003t ) 
found a similar behaviour in 4U 1728-34 ("Q2" in Fig. 9). 
Although it's premature to jump to conclusions with such 
a limited sample, these are at least good indications that 
NSXRBs could mimic the behaviour of some BHXRBs from 
this point of view. 

iMaitra fc Bailvnl (|2008t ) found a correlated behaviour 
between the optical (R band) and X-ray bands depending on 
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the state of the system. Using the same optical data set, and 
PCA X-ray data instead of ASM, we confirmed their find- 
ings. Moreover, we showed that trends exist also between the 
optical (R band) and HEXTE X-rays ba nds, and between 
the P CA and HEXTE bands (Fig. fO). iMaitra fc Bailvnl 
lj2008h proved that during the outbursts the evolution of 
the optical/infrared colour magnitude diagrams of Aql X- 
1 is consistent with thermal heating of an irradiated outer 
accretion disc. They noted that the optical/infrared colour 
and flux do not seem to change suddenly during the X-ray 
state transitions. This can be interpreted as evidence against 
a non-thermal origin of the optical/infrared emission. We 
didn't find any significant correlation between the optical (R 
band) and radio emission (at cm wavelengths), and although 
based on a very limited number of radio observations, this 
adds further support to this interpretation. 

As a general remark, our data are roughly consis- 
tent with a truncated disc mo del of the X-ray states (e.g. 
iMcClintock fc Remillardll2006l for a review). In hard X-ray 
states a standard accretion disc is truncated relatively far 
from the compact object. In soft X-ray states the disc is 
advancing closer to the innermost stable circular orbit. In 
the soft states the radiation is well described by black- 
body thermal emission from the inner regions of the accre- 
tion disc. But in the hard states the nature of the emit- 
ting material filling the gap between the compact object 
and the truncated disc (g enerically called the "corona" ) 
is sti l l under debate (e.g. lEsin et al.l I2001I : iMarkoff et al.l 
I2OOII : iBlandford fc BegelmanI 120041 ') . although it is quite 
clear that synchrotron and Compton mechanisms are im- 
portant. The fact that the optical emission (R band) does 
not seem to change significantly in the two X-ray states 
is further circum s tantia l evidence for the interpretation of 
iMaitra fc Bailvr] (|2008l ) of its origin as irradiation of the 
outer accretion disc. 
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